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We perform a measurement of the t —>■ I'yvv {I = e, p) branching fractions for a minimum photon 
energy of 10 MeV in the r rest frame, using 431 fb“^ of e'''e~ collisions collected at the center-of- 
mass energy of the T(45') resonance with the BABAR detector at the PEP-II storage rings. We find 
B(r ^ pjuO) = (3.69±0.03±0.10)xl0■^ and B(t ejup) = (1.847±0.015±0.052) x 10■^ where 
the first quoted error is statistical, and the second is systematic. These results are substantially 
more precise than previous measurements. 

PACS numbers: 13.30.Ce, 13.35-r, 13.40.Em, 13.40.Ks, 14.60.Fg 


Leptonic r decays are generally well suited to inves¬ 
tigate the Lorentz structure of electroweak interactions 
in a model-independent way [l[. In particular, leptonic 
radiative decays r —>■ IjizD, where the charged lepton (1) 
is either an electron (e) or a muon (/r), have been studied 
for a long time because they are sensitive to the anoma¬ 
lous magnetic moment of the r lepton 0 . At tree level, 
these decays can proceed through three Feynman dia¬ 
grams depending on whether the photon is emitted by 
the incoming r, the outgoing charged lepton, or the in¬ 
termediate W boson, as shown in Fig.[T] The amplitude 
for the emission of the photon by the intermediate boson 
is suppressed by a factor {rrir/Mw)'^ with respect to a 
photon from the incoming/outgoing fermions and is thus 
negligible with respect to next-to-leading order (NLO) 
QED radiative corrections [^. Both branching fractions 
have been measured by the CLEO collaboration. CLEO 
obtained B{t —>• p'jizD) = (3.61 ± 0.16 ± 0.35) x 10“^, 
and B{t — >• e^vD) = (1.75 ± 0.06 ± 0.17) x 10“^ for 
a minimum photon energy of 10 MeV in the r rest 
frame [J. In addition, the OPAL collaboration finds 
B{t -t p'yvT?) = (3.0 ± 0.4 ± 0.5) x 10“^ for a minimum 
photon energy of 20 MeV in the r rest frame [^. 
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In the present work we perform a measurement of 
T —>■ branching fractions for a minimum photon 

energy of 10 MeV in the r rest frame. This analysis 
uses data recorded by the BaBAR detector at the PEP- 
II asymmetric-energy e+e" storage rings operated at the 
SLAG National Accelerator Laboratory. The data sam¬ 
ple consists of 431 fb“^ of e+e“ collisions recorded at at 
the center-of-mass energy (CM) -y/i = 10.58 GeV Q. The 
cross section for r-pair production is arr = 0.919 ± 0.003 
nb 0 corresponding to a data sample of about 400 x 10® 
T-pairs. A detailed description of the BaBAR detector is 
given elsewhere |^. Charged particle momenta are 
measured with a five-layer double-sided silicon vertex 
tracker and a 40-layer helium-isobutane drift chamber 
inside a 1.5 T superconducting solenoid magnet. An 
electromagnetic calorimeter (EMC) consisting of 6580 
CsI(Tl) crystals is used to measure electron and pho¬ 
ton energies; a ring-imaging Cherenkov detector is used 
to identify charged hadrons; the instrumented magnetic 
flux return (IFR) is used for muon identification. About 
half of the data were taken with the IFR embedded with 
resistive plate chambers, later partially replaced by lim¬ 
ited streamer tubes. 

For this analysis, a Monte Carlo (MC) simulation is 
used to estimate the signal efficiency and to optimize the 
selection algorithm. Simulated r-pair events are gener¬ 
ated using KK2f (T^ and r decays are simulated with 
Tauola [^. Final-state radiative effects for r decays in 
Tauola are simulated using Photos [T^ . A signal r-pair 
MC sample is generated where one of the r leptons de¬ 
cays to r —>• Iji'D, and the other decays according to 
known decay modes [l^ . For the signal sample we re¬ 
quire the minimum photon energy in the r rest frame 
to be A* > 10 MeV. The r —^ IjvD decays with 
A* jnin < 10 MeV are treated as background. A separate 
r-pair MC sample is generated requiring each r lepton to 
decay in a mode based on current experimental knowl¬ 
edge; we exclude signal events in the former sample to 
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FIG. 1: Standard Model Feynman diagrams for r —^ lyvv at tree level. 


obtain a r-pair background sample. Other MC simu¬ 
lated background samples include qq {uu, dd, ss, 

cc), and BB {B = i?+, i?°) events. The events 

are generated by KK2f qq events are generated using 
the JETSET generator while BB events are simu¬ 
lated with EVTGEN lla. The detector response is sim¬ 
ulated with GEANT4 [2. Background from two-photon 
and Bhabha events is estimated from data. 

The signature for r —>■ l^vv decays is a charged par¬ 
ticle (track), identified either as an electron or a muon, 
and an energy deposit (cluster) in the EMC not associ¬ 
ated with any track, the photon. Since r leptons decay 
mostly to a single charged particle, events with two well- 
reconstructed tracks and zero total charge are selected, 
where no track pair is consistent with being a photon 
conversion in the detector material. The transverse mo¬ 
mentum of each track is required to be px > 0.3 GeV/c, 
the cosine of the polar angle is required to be between 
—0.75 and 0.95 within the calorimeter acceptance range 
to ensure good particle identification. The total miss¬ 
ing transverse moment of the event is required to be 
PT.miss > 0.5 GeV/c. All clusters in the EMG with no 
associated tracks (neutral clusters) are required to have 
a minimum energy of 50 MeV. We also reject events with 
neutral clusters having E < 110 MeV if they are within 
25 cm of a track, where the distance is measured on the 
inner wall of the EMC. 

Each event is divided into hemispheres (signal and tag 
hemispheres) in the CM frame by a plane perpendicu¬ 
lar to the thrust axis, calculated using all reconstructed 
charged and neutral particles [l^. Eor every event, the 
magnitude of the thrust is required to be between 0.9 
and 0.995. The lower limit on the thrust magnitude 
rejects most qq events while the upper limit removes 
e+e” —>■ and Bhabha events. The signal hemi¬ 

sphere must contain one track and one neutral cluster. 
The tag hemisphere must contain one track, identified 
either as an electron, muon or pion, and possibly one 
additional neutral cluster or n7r° (n =1, 2). Each 7r° 
candidate is built up from a pair of neutral clusters with 
a di-photon invariant mass in the range [100,160] MeV. 
To further suppress di-muon and Bhabha events, we re¬ 
ject events where the leptons in the signal and tag hemi¬ 
spheres have the same flavor. Since there are at least 
three undetected neutrinos in the final state we require 


the total energy to be less than 9 GeV. In the signal 
hemisphere, we require that the distance {di^) between 
the track and the neutral cluster, measured on the inner 
wall of the EMC, to be less than 100 cm. 

Electrons are identified by applying an Error Cor¬ 
recting Output Code (ECOC) algorithm based on 
Bagged Decision Tree (BDT) classifiers using as in¬ 
put the ratio of the energy in the EMC to the magnitude 
of the momentum of the track (E/p), the ionization loss 
in the tracking system (dE/dx), and the shape of the 
shower in the electromagnetic calorimeter. 

Muon identification makes use of a BDT algorithm, us¬ 
ing as input the number of hits in the IFR, the number of 
interaction lengths traversed, and the energy deposition 
in the calorimeter. Since muons with momenta less than 
500MeV/c do not penetrate into the IFR, the BDT also 
uses information the energy loss dE/dx in the tracking 
system to maintain a very low tt — p misidentification 
probability with high selection efficiencies. The electron 
and muon identification efficiencies are 91% and 62%, re¬ 
spectively. The probability for a tt to be misidentified as 
an e is below 0.1%, while the probability to be misiden¬ 
tified as a /X is around 1% depending on momentum. 

After the preselection, both samples are dominated 
by background events. For the r —> p'yvv sample, 
the main background sources are initial-state radiation 
(ISR), T —>• TTTT^zx decays, e+e“ —>■ p~ events, and 
r —>■ TTi/ decays. For the t —>■ e'yvv sample, almost 
all background contribution is from r —>■ evD decays in 
which the electron radiates a photon in the magnetic field 
of the detector (bremsstrahlung). Further background 
suppression is obtained by placing requirements on the 
angle between the lepton and photon in the CM frame 
(cosOij). For r —^ pji/P we require cos9i^ > 0.99, while 
for T —^ ejiyp we require cos^;.^ > 0.97 (see Figs. |2| 
and 13). To reject background from r —>■ evv decays in 
the T —>■ eyzxi/ sample, we further impose a minimum 
value for the invariant mass of the lepton-photon pair 
Mp >0.14 GeV/(? for this channel. In addition to the 
aforementioned quantities, the selection criteria use the 
energy of the photon and di^. The selection criteria are 
optimized in order to give the smallest statistical and 
systematic uncertainty on the branching fractions. 

After optimization, for r —> we require cos Op > 

0.99, 0.10 < Ej < 2.5 GeV, 6 < dp < 30 cm, and 
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Mi^ < 0.25 GeV/c^. The requirement on Mi^ rejects 
backgrounds from non-signal r decays. For the r —>■ e^vD 
channel, we require cos Oij > 0.97, 0.22 < < 2.0 GeV, 

8 < dij < 65 cm in addition to Mi^ > 0.14 GeV/c^. 

The signal efficiencies, the fraction of background 
events, and the number of events selected in the data 
are given in Table HI 

TABLE I: Signal efficiencies e (%), expected fractional back¬ 
ground contribution /bkg = A7,kg/(Afsig + Abkg), where Nsig 
is the number of signal events and A7,kg is the number of 
background events, and number of observed events (Aobs) for 
the two decay modes after applying all selection criteria. All 
quoted uncertainties are statistical. 



T —>■ fl')PU 

T eypu 

e 

0.480 ± 0.010 

0.105 ±0.003 

/bkg 

0.102 ±0.002 

0.156 ±0.003 

Nobs 

15688 ± 125 

18149 ± 135 


The branching fraction is determined using 

„ .^obs(l ~ /bkg) 

El = — t; - r - 

where fVobs is the number of observed events, UrT is the 
cross section for r pair production, C is the total inte¬ 
grated luminosity, and the signal efficiency e is deter¬ 
mined from the MG sample. 

After applying all selection criteria, we find 

= (3.69 ± 0.03 ± 0.10) X 10-3 
B{T^e-fvv) = (1.847 ± 0.015 ± 0.052) X 10-2 

where the first error is statistical and the second is sys¬ 
tematic. The systematic uncertainties on signal efficiency 
and on the number of the expected background events af¬ 
fect the final result, and are summarized in Table In] The 
most important contributions to the total uncertainty are 
from the uncertainties on particle identification, and pho¬ 
ton detection efficiency. 

To estimate the uncertainty on photon detection effi¬ 
ciency, we rely on e“''e- —>■ events for the high 

energy region [E^ > 1 GeV), and photons from 7r° de¬ 
cays for the low energy region [E^ < 1 GeV). Using fully 
reconstructed e+e- —?> events, we find that the 

photon detection efficiency for data and MG samples are 
consistent within 1% for U-j, > 1 GeV. For photon ener¬ 
gies E^ < 1 GeV, we measure the ratio of the branching 
fractions for r —?> ttz/ and t ^ pv decays. The resulting 
uncertainty on the reconstruction efficiency is found to 
be below 3%. Taking into account the 1.1% uncertainty 
on the branching fractions, the resulting energy-averaged 
uncertainty on the single photon detection efficiency is 
1.8%. We use this value as the systematic uncertainty in 
the efficiency for r —>■ Iji/D. 

The uncertainties on particle identification efficiency 
are estimated using control samples, by measuring the 


deviation of the data and MG efficiencies for tracks with 
the same kinematic properties. The uncertainty on the 
efficiency of the electron identification is evaluated using 
a control sample consisting of radiative and non-radiative 
Bhabha events, while the uncertainty for muons is an 
e“'"e- —>■ /X+/X-7 control sample. The uncertainty on the 
probability of misidentifying the pion as a muon or elec¬ 
tron is evaluated using samples of r —?> irinri' decays. The 
corresponding systematic uncertainty on the efficiency 
for r —>■ is 1.5% for both channels. 

For the background estimation, we define control re¬ 
gions that are enhanced with background events. For 
T —>■ p'yvv, where the major background contribution 
is not peaking in cosO^j, we invert the cut on cos0^^. 
For cos9fj,-y < 0.8, the maximum expected signal rate 
is 3% of the corresponding background rate. The maxi¬ 
mum discrepancy between the MG sample prediction and 
the number of observed events is 8%, with an excess of 
events in the MC sample. We take this discrepancy as es¬ 
timate of the uncertainty on the background prediction. 
For T —^ eqz/P, where the major background contribu¬ 
tions have similar cosOe-y distributions as signal, we ap¬ 
ply a similar strategy after requiring the invariant mass 
Mi^ < 0.14 GeV/c^; in this case we take cos0e.y < 0.90. 
The maximum contamination of signal events in this re¬ 
gion is 10%, and the maximum discrepancy between the 
prediction and the number of observed events is 4% with 
an excess of data events. We take this value as an es¬ 
timate of the uncertainty on the background rate. The 
errors on the branching fractions due to the uncertainty 
on background estimates are 0.9% for r —>■ p.'yvv, and 
0.7% for T —>• eqi/ix, respectively (Table HB. Gross-checks 
of the background estimation are performed by consider¬ 
ing the number of events expected and observed in differ¬ 
ent sideband regions immediately neighboring the signal 
region for each decay mode and found to be compatible 
with the aforementioned systematic uncertainties. 

All other sources of uncertainty, including current 
knowledge of the r branching fractions (BF), to¬ 
tal number of r pairs, limited MC statistics, dependence 
on selection criteria, and track momentum resolution are 
found to be smaller than 1.0%. 

In conclusion, we have made a measurement of the 
branching fractions of the radiative leptonic r decays 
T —>■ e'yi'v and r —>■ p.'^vv, for a minimum photon 
energy of 10 MeV in the r rest frame, using the full 
dataset of e+e- collisions collected by BaBAR at the 
center-of-mass energy of the T’(4S') resonance. We find 
B{t -a p-ivv) = (3.69 ± 0.03 ± 0.10) x 10-^, and B{t -a 
e^vv) = (1.847 ± 0.015 ± 0.052) x IQ-^, where the first 
error is statistical and the second is systematic. These 
results are more precise by a factor of three compared 
to previous experimental measurements. Our results are 
in agreement with the Standard Model values at tree 
level, B{t -a wjvv) = 3.67 x 10-3, and B{t -a ejvi?) = 
1.84 X 10-2 [3|, and with current experimental bounds. 
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FIG. 2: Selection of the r —^ fiyi/i?: (a) distance between lepton and photon candidates on the inner EMC wall, (b) cosine of 
the angle between momenta of the lepton and photon candidates in the CM frame, (c) invariant mass of the lepton photon 
pair, and (d) photon candidate energy in the CM frame for radiative r decay into a muon after applying all selection criteria 
except the one on the plotted quantities. The selection criteria on the plotted quantities are highlighted by the vertical lines; 
we retain the regions indicated by the horizontal arrows. 
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TABLE II: Summary of systematic contributions (%) to the 
branching fraction from the different uncertainty sources for 
the two signal channels. The total systematic uncertainties 
are obtained summing in quadrature the various systematic 
uncertainties for each decay channel. 



r ^ fx'yi'i' 

r 

Photon efficiency 

1.8 

1.8 

Particle identification 

1.5 

1.5 

Background evaluation 

0.9 

0.7 

BE [1^ 

0.7 

0.7 

Luminosity and cross section 

0.6 

0.6 

MC statistics 

0.5 

0.6 

Selection criteria 

0.5 

0.5 

Trigger selection 

0.5 

0.6 

Track reconstruction 

0.3 

0.3 

Total 

2.8 

2.8 
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FIG. 3: Selection of the r —>■ sample: (a) distance between lepton and photon candidates on the inner EMC wall, (b) 

cosine of the angle between momenta of the lepton and photon candidates in the CM frame, (c) invariant mass of the lepton 
photon pair, and (d) photon candidate energy in the CM frame for radiative r decay into an electron after applying all selection 
criteria except the one on the plotted quantities. The selection criteria on the plotted quantities are highlighted by the vertical 
lines; we retain the regions indicated by the horizontal arrows. 
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